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Abstract 



Identical particle correlations at fixed multiplicity are considered in the pres- 
ence of chaotic and coherent fields. The multiplicity distribution, one-particle 
momentum density and two-particle correlation function are obtained based on 
the diagrammatic representation for cumulants in semi-inclusive events. Our 
formulation is applied to the analysis of the experimental data on the multi- 
plicity dependence of correlation functions reported by the UAl and the OPAL 



Tj- . Collaborations. 
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^ '■ 1 Introduction 

a: 

In high energy hadron-hadron collisions, Bose-Einstein correlations of the identical 
particles are considered as one of the possible measures for the space-time domain 
where identical particles are produced. One of the theoretical approaches to the 
Bose-Einstein correlations is made on the analogy of the quantum optics Q, where 
two types of sources, chaotic and coherent are introduced. A diagrammatical method, 
based on the Glauber-Lachs formula |] , has been proposed to find the higher order 
Bose-Einstein correlation (BEC) functions in the quantum optical (QO) approach. In 
Ref. ||, the generating functional (GF) for momentum densities is derived in the QO 
approach, and a diagrammatic representation for cumulants is proposed. 

Up to the present, identical particle correlations in fixed multiplicity events are 
investigated in the case of purely chaotic field. Two-particle correlations are analyzed 
in Ref. Q by using Monte Carlo methods. Multiplicity dependence of one-particle 
distributions is discussed in Ref. , and that of two or three-particle correlations are 
considered in Ref. ||. 

In Ref. 0, an outline of our formulation on the particle correlations at fixed 
multiplicity in the QO approach has been briefly reported. General features of multi- 
plicity distributions, one-particle distributions and two-particle correlations at fixed 
multiplicity have been also shown. In the present paper, identical particle correlations 
at fixed mutiplicity in the QO approach are considered in detail. The diagrammatic 
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representation for cumulants is used to obtain the formulas in semi-inclusive events 
on the analogous way to that in inclusive events ||. Furthermore, our formulas are 
applied to the analyses of the experimental data in pp collisions by the UA1 Collab- 
oration [R], and in e + e~ collisions by the OPAL Collaboration |J. 

At first, we consider the case when there are no correlations among produced 
particles in the final states. Then particles in the final states are given by a coherent 
state, 

|*> = exp[~/ \f(p)\^ + J /(p) a+(p)^] |0). (1) 
The n-particle momentum density in semi-inclusive events is defined by 

Pn(Pi,---,Pn) = -^—E 1 ---E n d = |(0|a(pi) ■ • ■ a(p n )\4>)\ 2 , 

<j me \ & p\ ■ d p n 

which is reduced to 

Pn(Pn---,Pn) = \f(Pl)\ 2 ---\f(Pn)\ 2 eM- J Ifip)^}- (2) 

In the QO approach, the function f(p) is devided into two parts; 

M 

fip) = Y, a Mp) + fc(p), (3) 

where 4>i(p) and f c (p) are amplitudes of the ith chaotic source and a coherent source, 
and a,i is a random complex number attached to the zth chaotic source. In addition, 
M is the number of independent chaotic sources |J, which is regarded to be infinite 
in the present paper. The n-particle momentum density in the QO approach |L0| is 
defined by, 

Pn(PU---,Pn) = (\f(pi)\ 2 ---\f(Pn)\ 2 eM- J \f(P)\ 2 ^\j • (4) 

In Eq. (|4|), parenthesis (F) a denotes an average of F over the random number with 
a Gaussian weight ||; 

w-n(^-/«p[-^v)F. ( 5) 

It should be noticed that the classical (pion) fields are randomized in our ap- 
proach. On the other hand, each mode of the light is randomized in the quantum 
optics 0. After the average is taken over the random number in Eq. (|5]), terms 
composed of a^a*" 1 in the function F vanish if I ^ m. 

The generating functional (GF) for momentum densities in semi-inclusive events 
is defined by the following equation, 

Z sm [Hp)} = £ A /•••/ PnijPU- ■ ■ ,Pn)h{pi) ■ ■ ■ hijPn)^- • • • (6) 
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From Eqs. @ and (||), the GF is written formally as, 



Z SIa [h(p)) = ( exp 



/ \f(p) | 2 (%)-1)^ 



(7) 



On the right hand side of Eq.(||), an additional constant Z sra [h{jp) = 0], which does 
not affect to the momentum densities, is added. Inversely, the n-particle momentum 
density in the semi-inclusive events is given from the GF as 



Pn{pi 



,Pn) = Ei ■ ■ ■ E n 



S n Z sm [h(p)] 



Shipx) ■ ■■5h{p r , 



h(p)=0 



From the sum rule between semi-inclusive and inclusive cross-sections [nj, the 
GF Z[h(p)] for inclusive events is connected to that for semi-inclusive events; 

d 3 p 



Z[h{p)] = Z sm [h{p) + !] = ( exp 



fip) VKpY 



E 



The n-particle inclusive momentum density is given by 

5 n Z[h(p)} 



Pin(Pl,---,Pn) =E 1 ---E n 



Sh(pi) ■ ■■5(p n ) 



h(p)=0 



The explicit form of the GF, Eq. (§), is shown in Ref. ||, and the higher order 
BEC functions in inclusive events are obtained from it. In the following section, we 
would like to show that we can obtain higher order momentum densities in semi- 
inclusive events by analogy with a derivation in inclusive events. 



2 Generating functional and cumulant 

In the followings, we slightly change the definition of the GF from Eq.(^) to 



Z sm [h(p)] = c ( exp 



/ 1 f( P ) r 



(9) 



where the exponential damping factor in Eq. (0) is replaced by a normalization con- 
stant cq- Then, the n-particle momentum density in the QO approach is given by 



PniPl, ■•■,P2) 



E X -- - En 



5 n Z[h(p)} 



h(p)=0 



5h(pi) ■ ■■5{p n ) 
= c (\f(pi)---f(p n )\\. (10) 

The GF G sm [h(p)} for cumulants is defined by the following equation, 

G sm [h(p)]=lnZ sm [h(p)], (11) 



and the nth order cumulant is given by 



(12) 

fc(p)=0 



'5/i(pi) • • -5/i(p. 

From Eqs. (]TD|) , (ITTp and (O), we have iteration relations for momentum densities, 
P1O1) = c gx(pi), 

Pn(Ph---,Pn) = 9l(Pl)pn~l(p2,---,Pn) 
n-2 

+ Yl Y;9i+l{PliPhi ■ ■ ■ ,Pji)Pn-i-l(Pj i+1 , ■ ■ -,Pjn-l) 
i=l 

+ c g n (p u - ■ ■ ,p n ). (13) 

The second summation on the right hand side of Eq. (0) indicates that all possible 
combinations of (ji, • ■ ■ ,ji) and (ji+i, ■ ■ ■ ,j n -i) are taken from (2, 3, • ■ ■ , n). Equation 



(0) shows that the n-particle momentum density p n (pi, • • • ,Pn) ( n — 1> 2, ■ ■ •) can be 
evaluated if the cumulant g n (pi, • • • ,Pi) (i = 1, 2, • • • ,n) is obtained. 

The semi-inclusive one-particle and two-particle cumulants are given from Eq. ( |i~2|) 
respectively as 

9i(Pi) = (f(pi))a = r(pi,pi) + c(p 1 ,p 1 ), 
92(pi,P2) = (/(Pi)/(P2)>. - </(pi)>.</(pa)>« 

= k(pi,P2)| 2 + 2Re[r(pi,p 2 )c(p 2 ,Pi)], 

where r(pi,p 2 ) is a correlation caused by the chaotic sources and c(pi,p 2 ) is a corre- 
lation by the coherent source. Those are given by 

M 

r(pi,Pz) = ^Ai0i(pi)0*(p 2 ), 

i=l 

c(pi,Pa) = /c(pi)/ c *(P2). (14) 

As the GF, Eq.([|), in semi-inclusive events is the same with the GF, Eq.(H), 
in inclusive events except for the normalization factor cq, the cumulants of semi- 
inclusive events are also calculated from the same diagrammatic presentation as those 
of inclusive events [^J. Diagrammatic representation for them up to the fourth order 
is shown in Fig.l. 

The nth order cumulant is made up of connected terms of the correlations of 
chaotic fields, and those Cjj of the coherent field. If n > 3, the nth order cumulant 
gnijPii " " ' j Pn) is simply expressed by using the n-gon. The nth order cumulant consists 
of two types of terms; one is made only of the correlations of the chaotic fields, and 
the other contains the correlations both of the chaotic fields and of the coherent field. 
One of the terms belonging to the former type is given by ri 2 r 23 • • •r( n _i) n r n i. It can 
be expressed by the circular permutation (12 • • -n) started from 1. Other terms can 
be made from (12 ■ ■ - n) by any permutation. Therefore, there are (n — 1)! different 
terms in the former type. If any one of the correlations of the chaotic fields, 
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belonging to the former type, is replaced by the correlation c^- of the coherent field, 
the terms of the latter type are made. Therefore, there are n! different terms in the 
latter type. It should be noted that the half of the terms in the nth order cumulant 
are complex conjugates to the other half. 



3 Basic formulas at fixed multiplicity 

In order to calculate momentum densities at fixed multiplicity, the fc-particle momen- 
tum density and kth order cumulant at n-particle events (k < n) are defined by the 
following equations, respectively 



1 



(n-k)\ 
1 

(n - k,)\ 



PniPl,- ■ • ,Pk,Pk+l,- ■ ■ ,Pn) 
9n(Pl,- ■ ■ ,Pk,Pk+ir ■ ■ ,Pn] 



d 3 p k +i _ _ _ d 3 p n 
Ek+i E n 
d 3 p k+1 d 3 p Ti 



• (15) 



The GF for multiplicity distribution P(n) is given from Eq. 
independent from momentum p\ 



Z sm (h) = c (exp 



Ek+i E n 
@ if function h(p) is 

(16) 



The multiplicity distribution is given from Eq. (|T6|) ; 



= 1 8"ZUh) 



n\ dh n 



h=0 



I rn 

nl 



(17) 



The normalization of the fc-particle momentum density at n-particle events is given 
by 



p [ n\Pli---iPk) 



d 3 p 1 d 3 p k 



nl 



Pin). 



E 1 E k (n-k)\ 
Then the normalized fc-particle momentum density in n-particle events is defined as 

(n-k)\pW(pi,---,p h ) , 10 , 



P ( n\Pl,--- ,Pk) 



111 



Pin) 



In general, the inclusive fc-particle momentum density is given from the semi- 
inclusive momentum densities as 

,(*) 



pL (pi,---,Pk) = J2Pn ] (P^-- 



n=k 



which satisfies 



(k), ^d 3 pi d 3 p k ^ 

Pin '{Pl,---,Pk}— — = 2^ 



nl 



E x 



Ei 



^k(n-k)l 



P(n). 
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It should be noted that if the kth order inclusive momentum densities is integrated 
over all of the momenta, it becomes the fcth order factorial moment. 

From Eqs. fllTf), (|15D and ([[7]), the multiplicity distribution, and particle densi- 
ties in semi-inclusive events up to the third order can be expressed by the following 
recurrence equations, 



P{n) 
,(1)/ 



1 A 



i=i 



(19) 
(20) 
(21) 



i=i 

n— 1 n 
Pn ] {Pl,P2) = ^29j l) {Pl)pn-j{P2) + J29f\pi,P2)P{n- j), 
3=1 3=2 
n-2 

Pn ] (Pl,P2,P3) = ^9j(Pl)Pn-j(P2,P3) 
3=1 
n-1 

+ Yl {df (PuP2)Pn ) -j (Ps) + Of {Pl,P3)pn ) - j fe) } 
3=2 
n 

+ Y,9f\puP2,P3)P(n-j), (22) 

J'=3 

where -P(O) = cq. On the other hand, cumulants at fixed multiplicity are obtained 
from Fig.l as 



(0) = A[ R) +Ai s \ 



91 
9^ 



n-2 

Af)+2A2 1 + EAf n ) _,_ 1 

3=1 



n = 2,3, 



(23) 



9i ] (Pi) = r{p h pi) + c(p 1 ,p 1 ), 
9n\pi) = RnipuPi) + 25'„_i(pi,pi) + J2 T j,n-j-i(Pi,Pi), n = 2,3,---, (24) 



n-2 



i=i 



^ 2) (Pl>P2) 



n-1 

= R j(PuP2)Rn-j(P2,Pl) + 2c{p 1 ,p 2 )R n -l(P2,Pl) 

3=1 

n-2 

+ 2 J2{ S j(Pl>P2) R n-j-l(P2,Pl) + #n-j-l(Pl,P2)£j(P2, Pi)} 
3=1 
n-3 n-j-2 

+ E E {^v(Pl>P2)-Rn-j-i-l(P2,Pl) +-R„-j-i-i(pi,P2)^(p 2 ,p 1 )}, 
3=1 1=1 

(25) 



n— 2 n—j—1 

9n ] {Pl,P2,P3) = {Rj{Pl,P2)Rl{p2,P3)Rn-j-l{P3,Pl) 
3=1 1=1 
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where 



+Rj(PuPs)Ri (P3 ,P2)Rn-j-l(p2,Pl)} 
n-2 

+ 2 H {c(pi,P 2 )Rj (P2,P3)Rn-j-l(P3, Pi) 
J'=l 

+c(pi,p 3 )-Rj(P3,P2)-Rn-j-l(p2,Pl) 

+C(p 2 , P3 ) #n- j - 1 (P3 , Pi ) #j (Pi , P2 ) } 
n-3 n-j-2 

+ {5 , i (pi,p 2 )^(P2,P3) J Rn-i-«-l(P3,Pl) 

+ Sj (pi , p 3 ) (p 3 , P2 ) i2n - j - 1 - 1 (P2 , Pi ) 
+R j (p 1 ,p 2 )Si(p 2 ,P3)Rn-j-l-l(P3,Pl) 
+R j (p 1 ,p 3 )Si(p 3 ,p 2 )R n -j-l-l(p2,Pl) 
+R j (p 1 ,p 2 )R„-j-l-l(p2,P3)Sl(p3,Pl) 

+Rj(Pl,P3)Rn-j-l-l(P3,P2)Sl(p 2 ,Pl)} 

n—4 n—j—3 n—j—l—2 

+ E E 2 { 7 jv(Pl>P2)-Rm(P2,P3)-Rn-j-i-m-l(P3,Pl) 
j'=l 2=1 m=l 

(pi , p 3 ) it! m (p 3 , p 2 ) j_j_ m _ i {p 2 , Pi ) 
+R m (p 1 ,p 2 )T ji i(p 2 ,p 3 )R n _ j _i_ m _ 1 (p 3 ,p 1 ) 

+R m (pi, P3)Tj,l(p3, P2)Rn-j-l-m-l(P2, Pi) 
+Rm(Pl,P2)Rn-j-l-m-l(P2,P3)T j ,l(p 3 ,Pl) 

+R m (p 1 ,p 3 )R n -j-i- m -i(p3,P2)T j: i(p 2 ,p 1 )}, (26) 



Ri(pi,p 2 ) = r(p u p 2 ), 

/d^ k 
r (Pi, k)R n _i(k,p 2 ) , n = 2,3, 

S (Pl,P2) = C(p!,p 2 ), 

/d^ k 
c(p 1 ,k)R n (k,p 2 ) , n=l,2, 
to 

TjM,P2) = [ I RjipiMHhMM^) — ^-, (27) 

J J LUi U3 2 



Af) = jR n (k,k/ k 

= Js n (k,k)- 
Aj? = jT^k)^. (28) 



d 3 k 



In the followings, variables are changed from (pil, Pit) to (yi, Pit), with rapidity 
y! = tanh -1 ^!^/^!). Both correlations r(pi,p 2 ) and c(pi,p 2 ) are assumed to be real 



and parametrized as, 

r(yi,Prr',V2,Pzr) = P\f p(yi, P12X2/2, P2t) I(&y, Ap 1T ), 
c(j/i,Pit;2/2,P2t) = (1 - p)\/p(yi, Pir)p(l/2, P2T), 



7T 7T 



p(2/i,Pit) = < > \/-^ ex P [-« 2/i- Pir]> 

yap 

/(Ay, Ap T ) = exp[-7 L (A 2 /) 2 - 7T (Ap T ) 2 ], (29) 

where p = r(pi,pi)/p(pi) is called the chaoticity parameter, Ay = y 2 — y\ and Ap T = 
P2T — Pit- Functions defined by Eqs. fl2"7|) , (|2S|) and (^Uj) are expressed as 

Pit, 2/2, P2t) = ^ exp[-Aj(y 2 + y 2 2 ) + 2C 7 -y 1 y 2 ] 

x exp[-E^(p? T + p 2 T ) + 2W j p 1T ■ p 2T ), 

(1 - p) < n > a 1/2 /3 



Sj{yi,PiT,y2,Pzr) 



^Aj + a/2(Uj + (3/2) 



x exp 
x exp 



P 2 



a 



+ 



«7l 



2 Aj + a/2 



vi 



7i,j(?/l,PlT,J/2,P2T) 



2 4 + 0/2, 

(1-p) <n >7T 3 / 2 « 1 / 2 /? 



P2T 



^j(A i + a/2){A J + a/2){U l + 0/2) (E^ + 0/2) 



x exp 
x exp 



/ a 
-hr + 



"7l 



2 4 + ^/2, 

'0 07r N 
,2 4 + 0/2, 



P^T 



a 



"7l 



2 4 + a/2 y 

/0 ^7r_ X 
V2 4 + 0/2, 



y\ 



■2 

P2T 



(30) 



where 



A, 



4+1 



a 

2 +7L ' 



4- 



7l 



A,- + Ai 



2 



+ 7t, 



7t 



17.- + t/i 



W\ = 7 t, 



7Lfi- 
Aj + Ai ; 

7T^j 

Uj + Ui 



p < n > 



a 



X / 2 



iV 



7T 3 / 2 ' 

p <n > a 1/2 /3 



-iV„- 



A^ + Ai^ + C/O 



(31) 
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4 Analyses of experimental data 



Our formulas are applied to the analyses of negatively charged particles, or like- 
sign particles. Recently, the preliminary data of identical two-particle correlations in 
semi-inclusive events in pp collisions at \fs = 900 GeV within the pseudo-rapidity 
interval from -3.0 to 3.0 are reported by the UAl Collaboration ||. As can be seen 
from Eq. fl2T|) , multiplicity distribution and one-particle densities are included in the 
formula of the two-particle density, at least the multiplicity distribution of the same 
data sample is required to analyze the two-particle correlation. The multiplicity 
distribution at y/s = 900 GeV is also reported by the UAl group [O. However, the 
data are taken within the pseudo-rapidity interval from -2.5 to 2.5. In the present 
analysis, those data are used to adjust the parameters included in our formulation. 

Parameter a is determined from the parametrization of one-particle rapidity dis- 
tribution of Landau's hydrodynamical model [^3| , and (3 is taken from the inclusive 



transverse momentum distribution. In the present analysis , we neglect the correla- 
tion in the transverse momentum space, in other words, 7t is taken to be zero, and 
the calculated value is compared with the data after integrated over the transverse 
momentum. Therefore, the parametrization of f3 does not affect the calculated results. 
Those values are taken as 

a = 0.25, /5 = 5.556, 7 T = 0. 

Multiplicity distribution is normalized to satisfy 

E P{n) = 1, (32) 

71=1 

where n max = 42 is the maximum multiplicity of the observed negatively charged 
particles. 

Other parameters are adjusted to fit the multiplicity distribution from n = 1 



to n = 35 in the following way. The chaoticity parameter p is changed from p = 
to p = 1.0 by the step 0.1, and other parameters (no) and 7^ are determined by the 
minimum chi-squared method. The best fit is given by 

p = 0.8, (n ) = 2.584, 7l = 1.394, 

with Xmin/NDF = 372.4/32. As can be seen from Eq. (p9|), if we keep the relations 
that a/7i=costant and /3/jt— const ant, we get the same minimum chi-squared value. 
The calculated multiplicity distribution is compared with the data in Fig. 2. For the 
sake of comparison, we also fit the data by the negative binomial distribution, which 
results in xLn/ NDF = 465.0/33. 

The normalized one-particle rapidity distribution at n-particle events is defined 

by 

P [ n\y) = J ^Wt^Pt, (33) 

and calculated results for n = 5, 10 and 20 are shown in Fig.3. The peak height 
increases and the width of the distribution becomes narrower, as the multiplicity n 
increases. 
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In Fig. 4, the normalized two-particle rapidity distribution given by 



~p {2) (Ay) = J J j ~P [ n ] (Vi, Pit, yi + Ay, P2T)d yi d 2 p 1T d 2 p 2T (34) 

are shown at n—5, 10 and 20. The peak of the distribution also becomes higher and 
it's width becomes narrower, as the multiplicity n increases. However the increasing 
rate is gentler than that of the one-particle density. 

The normalized two-particle correlation function C^~\Ay) at n-particle events 
is defined as 



C (2 ~ ] (Ay)- ^ ^ P ( n ) ^^^,yi + Ay,p 2 T)dy 1 d 2 p 1T d 2 p2T 1 ^ 
/ / / Pn ) (Vi, Pir)Pn } (yi + Ay, p 2T )dy 1 d 2 p 1T d 2 p 2T 



The calculated result on C^~\Ay) at n=5, 10 and 20 are shown on Fig. 5. The 
multiplicity dependence of C!?~\Ay) at Ay=0 and 1.5 are shown in Fig. 6, where the 
preliminary experimental data reported by the UAl Collaboration at Q — 0.1 GeV 
are also shown []. 

In e + e~ collisions, the OPAL Collaboration published the data on multiplicity 



distributions [fPl], and multiplicity dependence of two-particle Bose-Einstein correla- 
tions |§ at 91 GeV. However, using the parameters adjusted to the observed mul- 
tiplicity distribution, which is close to a Poisson distribution, calculated results on 
C^~\Ay = 0) are almost constant and does not show the n dependence. Next, the 
multiplicity dependence of the Bose-Einstein correlations at Q = GeV, which is 
estimated from the data with Q > 0.05 GeV, is directly analyzed by our formula. We 
can fairly well reproduce the n dependence of the data with n max = 27, a = 0.125, 
p = 0.55 and ji = 10.0, if the minimum values of our calculated results on C^ 2 ~^(0) + 1 
are renormalized to 1. The result is shown in Fig. 7. 



5 Summary and discussions 

The analytical formulas of multiplicity distribution and particle densities in semi- 
inclusive events are derived from the generating functional GF in the presence of the 
chaotic and coherent fields. Formulas are applied to the analysis of the multiplicity 
dependence of two-particle correlations among identical particles in pp Collisions by 
the UAl Collaboration |§ and in e + e~ collisions by the OPAL Collaboration |§. 
In the formula of two-particle correlation, multiplicity distribution and one-particle 
densities in semi-inclusive events are contained. Therefore, to fix the parameters in 
our formulas, at least the observed multiplicity distribution is necessary. 

In pp collisions, we adjusted the parameters using the multiplicity distribution 
taken from the different data sample from those of the two-particle correlation. Our 
calculated results with the constant chaoticity parameter well reproduce the gross 

1 The data on the two-particle correlation are given by the variable 4- momentum transfer squared; 
Q = y— (pi — P2) 2 GeV. Q = corresponds to Ay = 0. Therefore we compare our calculated results 
at Ay = with the data at the smallest Q value (Q = 0.1) reported by the UAl Collaboration. 
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features of the multiplicity dependence of the data , inspite of the values of C^(0) 
are smaller than the data at Q = 0.1 GeV about 20%. Some part of the deviation will 
be deduced to the fact that the parameters are determined by fitting the multiplicity 
distribution within a different pseudo-rapidity interval. 

In e + e~ collisions, we analyze the data of the two-particle correlation without 
fitting the multiplicity distribution. We can explain the multiplicity dependence of 
the two-particle correlation at Q — Qmm GeV observed in the experiment, using the 
constant chaoticity parameter. 

Calculated results on the normalized two-particle correlation in semi-inclusive 
events show that the peak of the distribution becomes lower as the multiplicity in- 
creases, even if the chaoticity parameter p is constant. This behavior is similar to 
the data of two-particle correlations in pp collisions by the UAl Collaboration, and 
in e + e~ collisions by the OPAL Collaboration. 

In this paper, we analyze the data with the same values of chaoticity parameter p 
and the correlation length ji in rapidity space. Our present analyses indicate that the 
coherent component is not negligible, in other words, the values of chaoticity param- 
eters are smaller than 1. One of the possible candidates for the coherent component 



is a contribution from the decay products of long lived resonances ||15|| . Another 
possibility to reduce the value of chaoticity parameter is contanination [16, ||. For 
example, about 20% of like-sign particles are not pions in the OPAL Collaboration ||. 

When the colliding energy of incident particles increases as in the forthcoming 
RHIC experiment, thousands of identical particles can be produced in an event. Then, 
the production domain of those particles can be analysed precisely event by event. 
In general, the values of parameters will change according to the multiplicity. If the 
fitted values of chaoticity parameter or correlation length change suddenly at some 
multiplicity, it will be a possible signature that a threshold of a new phenomenon will 
open at that multiplicity. 
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Figure Captions 

Fig.l Diagrammatic representation for cumulants in semi-inclusive events. 

a) Contribution from the chaotic field, r(pi,p2), is shown by the solid line with 
arrow orienting from point 1 to 2. That of the coherent field, c(pi,p2), is shown 
by dotted line with arrow from 1 to 2. 

b) Diagram for (^(Pi,^)- 

c) Diagram for gz{p\, P21 Pz) ■ All permutations of (2,3) should be taken for (i, j). 

d) Diagram for g^pi, ■ ■ ■ ,P4). Those of (2,3,4) should be taken for (i, j, k) 
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Fig. 2 Multiplicity distribution observed in pp collisions [13 is analyzed by our 
formula. Parameters are determined by the minimum chi-squared method: 
p = 0.8, < n >= 2.584 and 7 L = 1.394. 



Fig. 3 Normalized one-particle rapidity distributions at fixed multiplicity calculated 
with p = 0.8, < n >= 2.584 and j L = 1.394. 

Fig. 4 Normalized two-particle rapidity distributions at fixed multiplicity calculated 
with p = 0.8, < n >= 2.584 and 7 L = 1.394. 

Fig. 5 Normalized two-particle correlation functions at fixed multiplicity calculated 
with p = 0.8 and 7 L = 1.394. 

Fig. 6 Multiplicity dependence of normalized two-particle correlations. Full circles 
show the data at Q=0.1GeV in pp collisions f|. Open circles and open squares 
are calculated with p = 0.8 and jl = 1.394. 

Fig. 7 Multiplicity dependence of normalized two-particle correlations. Full circles 
indicate the values at Q = GeV, estimated from the data with Q > 0.05 GeV 
in e + e~ collisions |§. Open circles are obtained from our calculation at Ay = 0. 
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